Supercurrent and noise in point contact between two different 

superconductors 

S. K. Yip 

Institute of Physics, Academia Sinica, Nankang, Taipei 11529, Taiwan 



Abstract 

We show that, for a quantum point contact between two superconductors 
of different gap magnitudes, Andreev bound states do not exist for certain 
phase differences and gap ratios. Continuum states may dominate the trans- 
port, and the supercurrent noise is qualitatively different from the case of 
equal gaps. 

PACS numbers: 74.80.Fp, 74.50+r 
With nanotechnology an emerging subject, there has been much effort in understanding 
electrical transport in quantum point contacts [1]. A particular interesting sub- field is the 
study of contacts between two superconductors. Their transport properties such as current- 
phase (J — 0) relationships and current- voltage dependences have been studied theoretically 
[2-4] as well as measured experimentally [5,6]. There are further theoretical works on the 
important issue of current fluctuations (noise). [4,7-9] All these theoretical papers are for 
the case where the two superconducting gaps |A| are equal. In the physical pictures devel- 
oped in these papers, a particularly important role is played by the Andreev bound states 
that form between the two superconductors. For a short junction with a single conduc- 
tion channel of transmission probability D, there is a pair of bound states at ±e& where 

r -i 1/2 

€b = |A| 1 — Dsin (0/2) . [10,11] Here is the phase difference between the two super- 
conductors. In equilibrium, the net supercurrent / is carried entirely by this bound state, 
and at zero temperature, I is given simply by / = — while at a finite temperature T, 

1 



/ is reduced by a factor tanh [^j arising from the occupation numbers at ±e b . The re- 
sponse of these bound states also play an important role in determining the non-equilibrium 
electron transport such as IV characteristics and the junction behavior under microwave 
irradiation. [2-4] 

These Andreev bound states also play a major role in determining the supercurrent 
noise of the point contact. [4,7-9] An important source of the equilibrium supercurrent 
noise arises from the fact that these bound states can either be occupied or empty due to 
thermal fluctuations. [7,8] As a result, there are two strong peaks of the current noise at 
zero frequency and u = le^. 

In this paper we study the bound state, current-phase relationship, as well as noise of a 
quantum contact between two superconductors with unequal gaps. We shall show that, in 
some regions of parameter space, the situation is dramatically different from the case where 
the gaps are of equal magnitude. For certain gap ratios and phase difference, the Andreev 
bound states do not exist at all. Instead of having bound state carrying all the current, the 
continuum states with energy between the two gaps can carry a substantial, sometimes even 
the dominating, part of the supercurrent. Interestingly the total supercurrent, and hence the 
current-phase relationship, is not significantly changed (except the overall magnitude) even 
in the case of unequal gaps. However, when the bound states do not exist, the supercurrent 
noise is qualitatively different. In particular, there is a substantial reduction in noise level 
for all frequencies below twice the smaller gap. The zero frequency and le^ peaks are entirely 
absent. It should be possible to perform measurements such as STM with a Pb tip on Al to 
verify the predictions in the present work. 

We shall then consider two s-wave superconductors in contact through a barrier with 
transmission coefncent D. Without loss in generality we shall take the left (right) super- 
conductor with smaller (larger) gap magnitude |Aj| ( |A r |), and take the phase difference 
4> between and n. One can easily adapt the quasiclassical formalism [12] to the present 
(one conduction channel) case, and obtain the Green's functions gi- r {p,£) and gi, r {p,t) for 
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positive (p) and negative (p = —p) momenta and at a point near and to the left and right 
(subscript / and r) of the barrier. It is convenient to express the answers [12,13] in terms of 
h,r = gi,r(p) + 9iAp) and di,r = 9i,r(p) - 9i,r{p)- 0ne obtains d[ = d r = d with 

iD 

d= 2^ i9r,oo,gi,oo] I V (1) 

and 

h = ((2 - D)g l>00 + Dg rtOB ) / V (2) 

where V(e) = 1 + ^ (g ljOD - g r>00 ) 2 . Here g ljr>0O = -7r -^==gp are the quasiclassical Green's 
function at positions far away from the interface into the left (right) superconductor. The 
above equations hold for both the retarded (R) and advanced (A) Green's functions with 
e — > e ± i0 + and we have left out their superscripts temporarily for simplicity. Here + is an 
infinitesimal positive number. Explicitly, we have 

; • n e^3(A r - Aj) - i|A r ||Aj|sin(0)f 3 

d = mD W) (3) 

where Q(e) = yj\ A^| 2 - e 2 ^/ 1 A r | 2 - e 2 V(e), i.e. 

Q(e) = (1 - ^Vl^P-^Vl^l 2 -^ " y ( e2 " |A r ||A z |cos(0)) (4) 
Andreev bound state exists when 

Q(e) = (5) 

for some |e| < \A t \ < \A r \. The solutions to eq (5) always come in ±e& pairs and we shall 
focus on €b > 0. Eq (5) can be rewritten as a quadratic equation in e 2 (by putting one of the 
terms of (4) to the other side and taking the square). It can be shown that one of the roots 
is larger than |A^| 2 and thus must be rejected. Thus for a given phase difference, there is 
at most one pair (±e&) of bound states. The remaining root must be substituted back into 
eq (5) to check that it is indeed a solution. We found that as the ratio in gap magnitude 
increases, there is an increasing region near = where the bound states cease to exist. 



Examples are shown in Fig. 1 and 2. The value of where the bound states first appear is 
given by <p c = cos~ 1 (| A ; |/| A r |) ( e = |A4 is a solution to eq (5) at cos0 = |Az|/|A r |). This 
critical phase difference C is independent of D, and <p c — > n/2 when \Ai\ « |A r |. 

To understand this <f> c , it is convenient [11] to imagine two small normal regions between 
the two superconductors surrounding the barrier (thus a S;NJN r S r junction. Here S^ r are 
the two superconductors, I the barrier and Nj jr are the two normal regions.). Quasiparticles 
and holes in N^ r are transmitted and reflected by the barrier I while Andreev reflected if 
they incident on the superconductors. Recall that for a particle (hole) incident on a super- 
conductor, the Andreev reflected hole (particle) acquires a phase —cos -1 (j^j) m addition 
to that from the phase of the order parameter =Farg(A). Now one can easily understand C 
for the case with perfect transmission. In this case one can ignore I and identify Nj jr . Let 
us take, without loss in generality, r = and (pi = 0. One can readily see that one has 
an Andreev bound state at e = |A;| at = C consisting of a hole propagating to the right 
and a particle propagating to the left. In this case there is no phase shift at the Andreev 
reflection at the left superconductor (since e = |Aj|), while the two contributions at the 
reflection at the right superconductor cancels (<f> r — cos -1 (^-|) = ). 

Note that for the above case, the phase shift does not cancel for an incident particle 
(reflected hole) to the right superconductor. Since reflection at the barrier I can in general 
generate this amplitude, it may be thus surprising at first sight that the bound state at 
e = |A/| still exists when D^l. However, a careful analysis (along the lines of [11]) shows 
that the bound state still exists but with the amplitude of the right moving particle in N r 
vanishes identically. The amplitude of the particle incident from the left of the barrier (in 
N;) and the amplitude of the particle incident from the right of it (in N r ) are in such a ratio 
that the outgoing, right moving particle amplitude in N r still vanishes exactly. 

The supercurrent is given by 




(6) 



where we have used the short-hand d R A = d R —d A and /(e) is the Fermi function. Explicitly, 



4 



we find 

eD\ A r \ | A;|sin0 1 



tanh( — M 



+ eD|A r ||A,|sin0 /-|A.| j _ (1 - f VKI 2 - e^e 2 - |A*| 2 tanh(e/2T) 



vr/i i|A ( | "(I - f ) 2 (|A r | 2 - e 2 )(e 2 - |Aj|2) + (| )2[|A r ||A z |cos^ - e 2 l 2 ( j 



where we have defined a^ r = y |A^ r | 2 — e 2 . The first term, due to the bound states, is 
to be included only when these states exist. I b can also be found from — Examples 
of J;, were plotted in Fig 1 and 2. however, is not the only contribution to the current 
(as in the case for equal gaps). d R ~ A is also finite for |A/| < e < |A r | (though still vanishes 
for |e| > |A r | ). This (and thus the second term of eq (7)) represents the contribution from 
continuum states with energies in the above range incident from the left superconductor and 
Andreev reflected from the right one. These states carry a finite current when <fi ^ 0. For 
low transmission, the continuum states contribution can actually dominate. Two examples 
of the various contributions to the current are shown in Fig. 3. It turns out that the 
dependence of the total current on is not much different from the case with equal gaps, 
except the magnitude is modified. 

We finally turn to the current noise of our junction. We shall study the quantity [14] 

S a (w) = J dte tuJt (< I(t)I(0) > - < I > 2 ) (8) 

S a represents the spectral density for absorption by the junction of energy quanta u. S a can 
be written (c.f. [16]) in terms of the envelope functions [17] for the Green's function near 
the interface (see Ref [18]) as 

^ 2 - de 



where 



p f elf 

SM = 1 r\ -(l-f(e + u))f(e)K(e) (9) 



K{e) = i{Tr 4 [f 3 C^ A (e + uj)r 3 C^ A (e) + r 3 C R Z A (e + u)r 3 C R Z A {e) 

-r 3 C R Z A (e + co)r 3 C R - A (e) - %C R ~ + A {e + u)r 3 C R Z A {e)]} (10) 
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This formula is valid on either side of the junction, and we shall choose to do this on the 

left (and correspondingly leave out the subscripts / in C"s.) C++ and C can be expressed 

in terms of s and d given earlier via C±± = ±| + ^(s; ± d). C+- and C-+ are given by 
(c.f. [16], see also [18]) C< + - = |r* (i - ^, i0O ) /D and (7_ + = -|r (i + ig ltOB ) /V. Here r is 
the reflection coefficient for normal quasiparticles incident from the left of the barrier, and 
|rf — 1 — D. 

With these C"s, the expression for S a can be written as eq (9) with K = K± + K 2 where 

Kfau) = —4(1 — D) Im 

+D 2 { (jii, e+u iii, e + Hi,e + uM,e I A^l 2 ) + (I -> r)} (11) 

+D(2-D) {(/ii, e +>r, e +/i r , e+ ^,e) + {Hl, e +u>Hr,e + Hr,e+uHl,e) ( | A; 1 1 A r |cOS0) } 

and 
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Here we have introduced the short hands //j jr (e) 
Im 



Im 



(12) 

and /v(e) = 



. On the right hand sides of eq (11) and (12), the subscripts e or e + uj 



specify where the functions are to be evaluated and all e's are understood to be e + ir/. Here 
is a damping coefficient. [7,8] 

An example of the current-noise is as shown in Fig 4. At increasing gap ratios, e& moves 
towards the gap edge \Ai\. Correspondingly the peaks at uj = and uj = 2e b decreases. 
When the bound states no longer exist, these two peaks are entirely absent and the noise 
becomes low for all frequencies below twice the smaller gap. This happens in particular for 
small phase differences and large gap ratios. Note also that the contribution from continuum 
states is strongly suppressed when T « |A/|. 

In conclusion, we have studied the Andreev bound states and supercurrent for a quantum 
point contact between two superconductors with different gap magnitudes. We showed in 
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particular that the current noise can be qualitatively different from the case where the 
superconductors are identical. 
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FIGURES 
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FIG. 1. Bound state energy eb (in unit of |A/|, left axis), and the bound state contribution to 
the current (in unit of e\Ai\/h, right axis) as a function of (f>. D = 0.9, while |A r |/|Az|'s are given 
by 1.0 (full line), 1.25 (dot-dashed), 2.0 (dotted) and 10 (long-dashed). 
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FIG. 2. Bound state energy (in unit of |A/|, left axis), and the bound state contribution to 
the current (in unit of e|A;|/?i, right axis) as a function of <fi. |A r |/|A;| = 2.0, while D's are given 
by 0.9 (full line), 0.5 (dot-dashed), and 0.2 (long-dashed). 
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FIG. 3. Current-phase dependence for (a): D = 0.9, |A r |/|Aj| = 2.0. T = 0, and currents are 
in units of Long-dashed: dotted: continuum contribution, symbols: total. The current 

for | A r | = |A;| is also shown (full line) for comparison, (b): Same as (a) except D = 0.1. Note the 
difference in scale for /. 
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FIG. 4. S (a;) for D = 0.9, = 0.5tt, T/|A { | = 0.3, and damping coefficient r] = 0.02|A,|. 
|A r |/|A/| are given by, dot-dashed: 1.0, dotted: 2.0, long dashed 10.0. S a in unit of e ^ , w in 
unit of | Aj | . 
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